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. Abstract. Low temperature zero- bias conductance through two side-coupled quantum 
' dots is investigated using Wilson's numerical renormalization group technique. A low- 
, temperature phase diagram is computed. Near the particle-hole symmetric point localized 
^ • electrons form a spin-singlet associated with weak conductance. For weak inter-dot cou- 
D ' pling we find enhanced conductance due to the two-stage Kondo effect when two electrons 
occupy quantum dots . When quantum dots are populated with a single electron, the 
system enters the Kondo regime with enhanced conductance. Analytical expressions for 
' the width of the Kondo regime and the Kondo temperature in this regime are given. 

s 

I Key words: Quantum dots, Kondo effect, Two-stage Kondo effect, Fano resonance 

■ 

C ■ 

O ' 

Q ' 1. Introduction 

t-h ! The interplay between the Kondo effect (which involves coupling of the 
^ local moment to the conduction electrons) and the magnetic ordering of 
. the moments has been studied in the context of bulk materials such as 
[ heavy- fermion metals (Jones and Varma, 1987). Recent advances in nan- 
CN ■ otechnology have enabled studies of transport through single as well as 
, coupled quantum dots where Kondo physics and magnetic interactions play 
| an important role at low temperatures. A double-dot system represents 
"j*^ ■ the simplest possible generalization of a single-dot system which has been 
| extensively studied in the past. Recent experiments demonstrate that an 
£h extraordinary control over the physical properties, such as the intra-dot 
i-^ . coupling, can be achieved in multiple dot systems (Jeong et al., 2001; Craig 
C I et al, 2004; Holleitner et al., 2002; Chen et al., 2004). This enables direct 
experimental investigations of the competition between the Kondo effect 
and the exchange interaction between localized moments on the dots. One 
manifestation of this competition is a two stage Kondo effect (Hofstetter 
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and Schoeller, 2003). Experimentally, it manifests itself as a sharp drop in 
the conductance vs. gate voltage (van der Wiel et al., 2002). 

The Fano resonance is a characteristic of noninteracting electrons for 
which the shape of the Fano resonance can be analytically determined. 
However, the influence of interactions on the appearance of Fano resonances 
remains an open question. It has been recently observed in experiments on 
rings with embedded quantum dots (Kobayashi et al., 2002) and quan- 
tum wires with side-coupled dots (Kobayashi et al., 2004). The interplay 
between Fano and Kondo resonance was investigated using equation of 
motion (Bulka and Stefanski, 2001; Stefanski et al., 2004) and Slave boson 
techniques (Lara et al., 2004). 

We study a double quantum dot (DQD) in a side-coupled configuration 
(Fig. 1), connected to a single conduction-electron channel. Similar systems 
were studied previously with different techniques, such as: non-crossing ap- 
proximation (Kim and Hershfield, 2001), embedding technique (Apel et al., 
2004), and slave-boson mean field theory (Kang et al., 2001; Lara et al., 

2004) . Numerical renormalization group (NRG) calculations were also per- 
formed recently (Cornaglia and Grempel, 2005), where only narrow regimes 
of enhanced conductance were found at low temperatures and no Fano res- 
onances were reported. We will present a low-temperature phase diagram 
of a DQD as a function of intra-dot coupling strengths and gate-voltage 
potential indicating regions with enhanced conductance due to Kondo effect 
and regions of nearly zero conductance. In particular, when the intra-dot 
overlap is large and DQD occupancy is one, wide regimes of enhanced 
conductance as a function of gate-voltage exist at low temperatures due 
to Kondo effect, separated by the regimes where localized spins on DQD 
are antiferromagnetically (AFM) coupled. Kondo temperatures Tjc follow a 
prediction based on the poor-man's scaling and Schrieffer- Wolff transforma- 
tion. In the limit when the dot a is only weakly coupled, the system enters 
the "two stage" Kondo regime (Vojta et al., 2002a; Cornaglia and Grempel, 

2005) , where we again find a wide regime of enhanced conductivity under 
the condition that the high- and the low- Kondo temperatures (Tk and T^, 
respectively) are well separated and the temperature of the system T is in 
the interval Tg- < T «% 




Figure 1. Side-coupled configuration of quantum dots. 
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2. Model and Method 

The Hamiltonian that we study reads 

H = 5(n d -l)+5(n a -l)-t d J2(dta a + ald a ) 

a 

+ |(n d -l) 2 + ^(n a -l) 2 (1) 

+ J2 e fc c L c *' + V d( k ) ( c la d * + d a c ka) 
ka ka 

~~ Jad^a " &di 

where rid = J2 a d\d a an d n a = J2a a t a a- Operators d\ and aj are creation 
operators for an electron with spin a on site d or a. On-site energies of the 
dots are defined by e = 5—U/2. For simplicity, we choose the on-site energies 
and Coulomb interactions to be equal on both dots. Coupling between 
the dots is described by the inter-dot tunnel coupling t d . Dot d couples 
to both leads with equal hopping t' . Operator c ka creates a conduction 
band electron with momentum k, spin a and energy e& = — D cos k, where 
D = It is the half-bandwidth. Spin operator S = J2s s' c \®s,s' c s' is defined 
using Pauli matrices and J ad represents additional Heisenberg exchange 
interaction. The momentum-dependent hybridization function is Vd(k) = 
— (2/y/N + 1) t' sink, where N in the normalization factor is the number of 
conduction band states. 

We use Meir-Wingreen's formula for conductance in the case of pro- 
portionate coupling (Meir and Wingreen, 1992), which is known to apply 
under very general conditions (for example, the system need not be in 
a Fermi-liquid ground state), with spectral functions obtained using the 
NRG technique (Wilson, 1975; Krishna-Murthy et al., 1980; Costi et al., 
1994; Costi, 2001; Hofstetter, 2000). At zero temperature, the conductance 
is 

G = G 7rT Pd (0), (2) 

where Go = 2e 2 /h, pd(oj) is the local density of states of electrons on site 
d and Y/D = (t'/t) 2 . 

The NRG technique consists of logarithmic discretization of the con- 
duction band, mapping onto a one-dimensional chain, and iterative diag- 
onalization of the resulting Hamiltonian (Wilson, 1975). Only low-energy 
part of the spectrum is kept after each iteration step; in our calculations 
we kept 1200 states, not counting spin degeneracies, using discretization 
parameter A = 1.5. 
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3. Strong intra-dot coupling 



Figure 2. Conductance 
and correlation functions 
of DQD vs. 5. Besides dif- 
ferent values of td, indi- 
cated in the figure, other 
parameters of the model 
are T/D = 0.03, U/D = 1 
and J a d = 0. Temperature 
T is chosen to be far below 
T K , i.e. T<T K . 

In Fig. 2a we present conductance through a double quantum dot at 
different values of intra-dot couplings vs. 5. Due to formation of Kondo 
correlations, conductance is enhanced, reaching the unitary limit in a wide 
range of 5. Regimes of enhanced conductance appear in the intervals ap- 
proximately given by 5± < \5\ < 62, where 5i = td(2^/l + (U/4td) 2 — 1) and 
62 = (U/2 + td)- These estimates were obtained from the lowest energies of 
states with one and two electrons on the isolated double quantum dot, i.e. 
E 1 = U/2-S- t d and E 2 = -2t dy /l + (U/U d ) 2 + U/2, respectively. 

Using the above estimates we present a phase diagram in Fig. 3. In the 
gray region, called Kondo regime, with border lines given by 5i(t d ), ^{td), 
the Kondo effect is responsible for an enhanced conductance. Kondo plateaus 
in Fig. 2 fall in this regime. The Conductance is zero everywhere outside 
this region, except in the limit when t d — * 0, where a two-stage Kondo effect 
two stage is responsible for enhanced conductance, as discussed further in 
the text. 

To gain further physical insight, we focus on various correlation func- 
tions, defined within the DQD system. In Fig. 2b we show S, calculated 
from expectation value (S% ot ) = S(S + 1), where S to t = 8,1 + 8^ is the total 
spin operator. 5 reaches value 1/2 in the Kondo regime where G/Gq = 1. 
Enhanced conductance is thus followed by the local moment formation. 
This is further supported by the average double-dot occupancy (n) , where 
n = n a + nd, which in the regime of enhanced conductivity approaches odd- 
integer values, i.e. (n) = 1 and 3 (see Fig. 2c and 3). Transitions between 
regimes of nearly integer occupancies are rather sharp; they are visible as 
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Figure 3. Phase diagram of a DQD for U/D — 1 and J ac i = 0, obtained using analytical 
estimates as given in the text. Grey areas represent Kondo regimes where S ~ 1/2, 
(n) ~ 1, and G/Go ~ 1. In shaded area, called spin - singlet regime, where S ~ 
and (n) ~ 2, there is enhanced spin-spin correlation function, i.e. (S a ■ So) ~ —0.5 and 
G/Go ~ 0. The two-stage Kondo regime is explained further in the text. 



regions of enhanced charge fluctuations measured by Are 2 = (re 2 ) — (n) 2 , as 
shown in Fig. 2d. Finally, we show in Fig. 2e spin-spin correlation function 
(S a • Srf). Its value is negative between two separated Kondo regimes where 
conductance approaches zero, i.e. for —Si < 5 < 5%, otherwise it is zero. 
This regime further coincides with (n) ~ 2. Each dot thus becomes nearly 
singly occupied and spins on the two dots form a local singlet (S=0) due 
to effective exchange coupling J e g = At^/U. 




Figure 4- Kondo 
temperatures Tk vs - S as 
measured from the widths 
of Kondo peaks obtained 
from NRG calculations 
(full circles and squares). 
Analytical estimate, 
Eq. 3, is shown using 
dashed lines. The rest of 
parameters are identical to 
those in Fig. 2. 



6 



J. BONCA 



In Fig. 4 we present Kondo temperatures Tk vs. 5 extracted from the 
widths of Kondo peaks. Numerical results in the regime where (n) ~ 1 and 
3 fit the analytical expression obtained using the Schrieffer-Wolf transfor- 
mation that projects out states with even electron occupancy and leads to 
an effective single quantum dot problem with renormalized parameters. We 
obtain 

T K = 0.182^7 \/~PqJ exp[—l/p J] (3) 
with pqJ = ^{a/\Ei\ + (3/\E 2 - Ei\), where a = 1/2 and 



[i&* + u(y + Ji&* + u*)) 



The prefactor 0.182^7 in Eq. 3 is the effective bandwidth. The same effective 
bandwidth was used to obtain Tk of the Anderson model in the regime 
U < D (Krishna-Murthy et al., 1980; Haldane, 1978). 



4. Weak intra-dot coupling 



We now turn to the limit when td — ► 0. Unless otherwise specified, we choose 
the effective temperature T to be finite, i.e. T ~ 10~ 9 D, since calculations 
at much lower temperatures would be experimentally irrelevant. In this 
case one naively expects to obtain essentially identical conductance as in 
the single-dot case. As 5 decreases below <5 ~ U/2, G/Go indeed follows 
result obtained for the single-dot case as shown in Fig. 5a. In the case of 
DQD, however, a sharp Fano resonance appears at <5 = U /2. This resonance 
coincides with the sudden jump in S, (n), as well as with the spike in 
An 2 , as shown in Figs 5b,c, and d, respectively. The Fano resonance is a 
consequence of a sudden charging of the nearly decoupled dot a, as its level 
e crosses the chemical potential of the leads, i.e at e = 0. Meanwhile, the 
electron density on the dot d remains a smooth function of 5, as seen from 
(rid) i n Fig 5c. With increasing td, the width of the resonance increases and 
at td ^ 0.1, the resonance merges with the Kondo plateau and disappears 
(see Fig. 2a). The regime where the Fano resonance exists is also specified 
in the phase diagram, Fig. 3. 

We now return to the description of the results presented in Fig. 5a 
in the regime where 5 < U/2. As 5 further decreases, the system enters 
a regime of the two-stage Kondo effect (Cornaglia and Grempel, 2005). 
This region is defined by J c ff < Tk (see also Fig. 5f), where Tk is the 
Kondo temperature, approximately given by the single quantum dot Kondo 
temperature, Eq. 3 with p J = ^(1/|<5 - U/2\ + 1/|<5 + U/2\). This regime 
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Figure 5. Conductance and correlation functions at td/D — 0.001 (a,...,e). Dashed 
lines represent in a) G/Go and in b) S of a single quantum dot with otherwise identical 
parameters. Dashed line in c) represents (rid) of DQD. In f) a schematic plot of different 
temperatures and interactions is presented as explained in the text. NRG values of the 
gap in pd(<jj) at ui — and T <C are presented with open circles. Values of J e s and 
analytical results of T% are presented with dashed lines. For analytical estimates of T% 
a — 2.2 was used. Other parameters of the model are T/D = 0.03, U /D — 1 and J a d = 0. 



is also indicated in the phase diagram, presented Fig. 4, where condition 
J c fj = Tk was used to separate two-stage Kondo from the spin-singlet 
regime. Just below 5 < U/2, T falls in the interval, given by T K <C T <C Tk, 
where T K ~ Tk exp(— aTx/Jes) denotes the lower Kondo temperature, 
corresponding to the gap in the spectral density pdioj) at to = and a is 
of the order of one (Cornaglia and Grempel, 2005). Note that NRG values 
of the gap in pd(w) (open circles), calculated at T <C T K , follow analytical 
results for T K (S) when J e fj < Tr-, see Fig. 5f, while in the opposite regime, 
i.e. for J c ff > Tk, they approach J e ff. 

As shown in Fig. 5a for 0.3D < 5 < U/2, G/G calculated at T = W- 9 D 
follows results obtained in the single quantum dot case and approaches 
value 1. The spin quantum number S in Fig. 5b reaches the value S ~ 0.8, 
consistent with the result obtained for a system of two decoupled spin- 1/2 
particles, where (S 2 ) = 3/2. This result is also in agreement with (n) ~ 2 
and the small value of the spin-spin correlation function (S • S^), presented 
in Fig. 5c and 5e respectively. 

With further decreasing of 5, G/Gq suddenly drops to zero at 5 ^ 0.3D. 
This sudden drop is approximately given by T ^ T K (5), see Figs. 5a and f. 
At this point the Kondo hole opens in Pd{u) at u = 0, which in turn leads 
to a drop in the conductivity. The position of this sudden drop in terms of 
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5 is rather insensitive to the chosen T, as apparent from Fig. 5f. 

Below 5 S= 0.25D, which corresponds to the condition J c g ~ Tk(5), also 
presented in Fig. 5f, the system crosses over from the two stage Kondo 
regime to a regime where spins on DQD form a singlet. In this case S 
decreases and (S a Sd) shows strong anti-ferromagnetic correlations, Figs. 5b 
and e. The lowest energy scale in the system is J e fr, which is supported by 
the observation that the size of the gap in p(to) (open circles in Figs. 5f) is 
approximately given by J e g. 



5. Ferromagnetic coupling 



In Fig 6 we present results of conductance and correlation functions for 
different values of the ferromagnetic exchange coupling J ac [, td/D = 0.3 and 
U/D = 1. At finite td and J a d = the ground state of an isolated DQD 

containing two electrons is a spin singlet state. At J^/D = 2{^J \t\ + U 2 — 
U) ~ 0.33 spin singlet and triplet states become degenerate. With further 
increasing J a d > J^ d we expect formation of a S = 1 state in the regime 
when (n) ~ 2. Since in this case 2S > K where K = 1 is the number 
of channels, this systems falls into a class where the spin on the DQD is 
not fully compensated by the conduction electrons (Vojta et al., 2002b; 
Hofstetter and Schoeller, 2003). Results at finite J a d, presented in Fig 6, 
can be divided into two groups. For J a d < J^ d the main effect of increasing 
J a d is seen as expansion of the S = 1/2 Kondo regime where (n) ~ 1. For 
Jad > Jad-> * wo Kondo plateaus appear, one associated with 5 = 1/2 Kondo 
regime and the other with S = 1 Kondo regime. In both cases conductance 
reaches unitary limit, while in the transition regime it drops to zero. 




Figure 6. Conductance 
and correlation functions 
of DQD vs. 8 for various 
values of J a d and td/D = 
0.3. Other parameters of 
the model are identical to 
those in Fig. 2. Tempera- 
ture T is chosen to be far 
below T K , i.e. T < T K . 
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6. Conclusions 

In this work we have explored different regimes of the side-coupled DQD. 
Conclusions can be summarized as follows: a) when quantum dots are 
strongly coupled, wide regions of enhanced, nearly unitary conductance 
exist due to the underlying Kondo physics. Analytical estimates for their 
positions, widths, as well as for the corresponding Kondo temperatures, are 
given and numerically verified. When two electrons occupy DQD, conduc- 
tance is zero due to formation of the spin-singlet state which is effectively 
decoupled from the leads, b) In the limit when quantum dots are weakly 
coupled the Fano resonance appears in the valence fluctuation regime. Its 
width is enhanced as a consequence of interactions which should facilitate 
experimental observation. Unitary conductance exists when two electrons 
occupy DQD due to a two-stage Kondo effect as long as the temperature 
of the system is well below Tk and above T^. The experimental signature 
of the two-stage Kondo effect in weakly coupled regime should materialize 
through the inter-dot coupling sensitive width of the enhanced conductance 
vs. gate voltage. 

When intra-dot ferromagnetic coupling exceeds a critical value, there 
is a phase transition from a spin-singlet, non-conducting regime to S = 1 
Kondo regime where the spin on the DQD is under-screened but neverthe- 
less the conductivity reaches the unitary limit. 
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